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ABSTRACT 

We explore the ability of high energy observations to constrain orbital pa- 
rameters of long period massive binary systems by means of an inverse Comp- 
ton model acting in colliding wind environments. This is particular relevant for 
(very) long period binaries where orbital parameters are often poorly known from 
conventional methods, as is the case e.g. for the Wolf-Rayct (WR) star binary 
system WR 147 where INTEGRAL and MAGIC upper limits on the high-energy 
emission have recently been presented. We conduct a parameter study of the 
set of free quantities describing the yet vaguely constrained geometry and re- 
spective effects on the non-thermal high-energy radiation from WR 147. The 
results are confronted with the recently obtained high-energy observations and 
with sensitivities of contemporaneous high-energy instruments like Fermi-LAT. 
For binaries with sufficient long periods, like WR 147, gamma-ray attenuation 
is unlikely to cause any distinctive features in the high-energy spectrum. This 
leaves the anisotropic inverse Compton scattering as the only process that re- 
acts sensitively on the line-of-sight angle with respect to the orbital plane, and 
therefore allows the deduction of system parameters even from observations not 
covering a substantial part of the orbit. 

Provided that particle acceleration acts sufficiently effectively to allow the 
production of GeV photons through inverse Compton scattering, our analysis 
indicates a preference for WR 147 to possess a large inclination angle. Otherwise, 
for low inclination angles, electron acceleration is constrained to be less efficient 
as anticipated here. 

Subject headings: Stars: early-type - Stars: binaries - Stars: winds, outflows - 
Gamma rays: theory - Radiation mechanisms: non-thermal - WR 147 
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Introduction 



Massive early type stars (0-, early B-, Wolf-Rayet (WR) stars) are hot (Tcs > 10000 K), 
and possess some of the strongest sustained winds among galactic objects. Among the 
massive early-type stars, WR-stars have the highest mass loss rate M ~ 10~^'"~^ Mq/jt. 
Winds of massive stars a re generally sup ersonic, reaching terminal velocities of typically 
Voo > 1000 — 5000 km/s (ICassinelli 19791 ). Their kinetic energy, = Mt>^/2, however, 
rarely exceed 1% of the bolometric radiative energy output of typically ~ 10'^^ erg/s. When 
two such stars form a system, their winds collide, creating an interaction region of high 
temperature (T ~ 10'' "^K) shocked gas. The position of this collision region is determined 
by ram pressur e balance, and two (re verse) shocks follow where Fermi acceleration has been 



suggested (e.g.. 



Eichler fc Usov 1993h . 



Convincing observational evidence for particle acceleration to relativistic energies 
medi ated by their supers onic winds comes from the detection of non-thermal radio emission 



[e.g., 



Abbott et al. 19861 ). This component has been interpreted as synchrotron emission 



on the basis of the measured spectra (softer than the canonical value ~ 0.6, F^, oc z/"'') 
and high brightness temperatures of ~ 10^^''K, and serves also as a proof for the existence 
of magnetic fields in the shock environment. These findings triggered the development of 
theoretical models of non-thermal photon production in the shock region which is exposed 
to a strong radiation field in the UV range and high wind particle densities from the 
participating hot stars, and their magnetic fields. In these scenarios non-thermal high-energy 
emiss ion from the sho c k region, extending well into the 7-ray domain, has been predicted 



Montmerle 1979 



Eichler fc Usov 1993 



Romero et al. 1999 



Benaglia fc Romero 2003 



Miicke fc Pohl 2002 



Reimer et al. 2006 



Pittard fc Doughertv 20061 ). 



Generally, 7-rays are produced either by inelastic nucleon-nucleon interactions if 
the cosmic ray nucleons reach sufficient high energies to overcome the threshold of this 
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process, or are produced through leptonic processes hke inverse Compton scattering and 
non-thermal bremsstrahlung. Typically, inverse Compton produced fluxes are found to 
dominate over bremsstrahlung produced ones here. For close binary systems it has also been 
proposed that nucleons are accelerated in the shock up to EeV e nergies, with su bsequent 



photodissociation reactions and/or photomeson production (e.g.. 



Bednarek 20051 ) that may 



initiate pair cascades, eventually leading to detectable 7-ray flux levels. 

In the course of the last decade, several high energy phenomena have been 
observed supporting the occurrence of 7-ray production in massive star systems: 
the EGRET detections of 7-ray sources that are, although presently unidenti- 
fied, found in intriguing positional coinc idence with massive star binary systems 



f lKaul fc Mitra 1997 



Romero et al. 19991 ). binary systems containing a massive star 



and a compact companion have been discovered as sources of TeV-emission (e.g., 



Aharonian et al. 2005a 



Aharonian et al. 2005b 



Albert et al. 2006 



Albert et al. 20071 ) 



the southern massive star cluster Westerlund 2 shows extended TeV-emission 



( lAharonian et al. 2007i ). as does Cyg 0B2 (lAharonian et al. 2005c 



Albert et al. 20081 ). 



with winds of its massive stars likely involved in powering the emission. Recently, the 
very-long period bi nary system WR 1 47, consisting of a WN8(h) and a BO. 5 V star. 



and with ~ 650pc ( iMorris et al. 2000l ) the second closest of all Wolf-Rayet binaries, 
has been observed with the MAGIC telescope for 30.3 hours of good data. From these 
observations the MAGIC collaboration reported upper limits of -F>80GeV < 1.1 • 10~^^ cm~ ^ 
s~\ i^>330GeV < 3.1 ■ 10"^^ cm'^ s"^ and F>6oogcV < 7.3 • 10^^^ cm-^ s-^ ljAliu et al. 20081 ). 



Because of WR 147's proximity this system belongs to the few where high-resolution 
radio observations lead to resolving the collision region between the main sequence stars 
as an extended, slight 



( iDougherty et al. 1996 



y elongated non-thermal feature on VLA and MERLIN images 



Contreras et al. 19971 : 



Dougherty et al. 20001 ) in addition to the 



free-free emission from the spherical wind of the stars. Furthermore, X-ray imaging 
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with Chandra revealed an extended region at the loc ation of ram pressu re balance 



and was interpreted as the wind-wind collision zone (jPittard et al. 20021 ). However 



despite a multitude of observations dedicated towards this object, as well as its relatively 
short distance to Earth, the geometry of this system is not well known. The unknow ns 
start with WR 147's period, being at least 1350 years (ISetia Gunawan et al. 2001al ). 
the lack of clues about its argument of periastron, and the yet pending observational 
hints towards a determination of its orbital eccentricity. Proposed inclination angles, 
mainly deduced from mod eling of the extended radio and X-ray ima ges, cover a wide 



range: from "low" (< 30°: 



WiUiams et al. 1997 



Doughertv et al. 2003 



Contreras et al. 19991: 



Pittard et a. 



. 2006h to "large" (> 40°: 



Pittard et al. 20021 ) with respect to the plane 



of the sky. Consequently, any predictions for the 7-ray emission from this source must be 
concurrently evaluated in conjunction with the assumed parameter space when confronting 
with observational high energy measurements. 

The aim of this work is to scientifically explore and relate the sensitivity of 7-ray 
detections or limits by operating (H.E.S.S. /VERITAS, MAGIC, AGILE, Fermi Gamma-ray 
Space Telescope) instruments to the geometry (inclination, eccentricity, orbital phase) 
of the massive binary system WR 147. For this purpose we are using the emission 
model for colliding winds of long-period massive stars of Reimer et al. (2006), which 
self-consistently calculates the non-thermal particle population and its radiation after 
diffusive acceleration of particles taken from the pool of thermal wind material. This 
model includes all relevant radiative processes/losses, expansion losses and propagation 
effects such as particle convection out of the wind's collision zone for a generally not too 
optically thick environment, which is shown to be appropriate for WR 147. As indicated 
in Reimer et al. (2006) inverse Compton (IC) emission well dominates over all competing 
7-ray production mechanisms in wide massive binary systems, thereby determining their 
high-energy appearance. Hence we restrict our calculations to this process only. Both 
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Klein-Nishina effects as well as the angle-dependence of the cross section are taken into 
account, with the latter playing a decisive role for the results of our exploration. 

In Sect. 2 of this paper we provide an overview of the system parameters and 
measurements of the non-thermal emission component of WR 147. A brief reminder about 
the main characteristics of the non-thermal emission model of Reimer et al. (2006) for 
the collision region of long-period massive binary star systems is given in Sect. 3. Sect. 4 
is devoted to a parameter study of the set of free geometry describing quantities and 
their effects on the predicted non-thermal high-energy radiation from WR 147. These 
are confronted with all information gathered so far from measurements of WR 147's 
non-thermal emission component. We conclude in Sect. 5 with our main results of this 
work, the constraint on the geometry of the orbit on WR 147. 



2. The knowns (and unknowns) of WR 147 



Located somewhat off the galactic plane and in spatial coincidence with the Cygnus 
0B2 region, although n ot with the nearby unidentified TeV source TeV 2032+41 



(lAharonian et al. 2005d ). WR 147 is one of the radio-brightest massive star systems known. 



and, after 7^ V elorum, the second clos est WR-star to Eart 



1. Distance estima tes range from 



(630 ± 70) pc fIChurchwell et al. 19921 ) to (650 ± 130) pc flMorris et al. 2000h . We will use 



650 pc as WR 147's distan ce in the following. A lthough the binary nature of this system was 



already proposed in 1990 (IWilliams et al. 19901 ) . it was not generally accepted until high 



resolution images undoubtedly established the existence of a weak Wolf-Rayet companion of 
type WN8(h) to the main sequence star BO. 5 V. In particular, WR 14 7 is meanwhile resolve d 



into it s components in the rad i o band by MERLIN 



VLA fIChurchwell et al. 1992 



Skinner et al. 199 



(1.6 GHz, 5 GHz) flWiUiams et al. 19971 ). 



Contreras et al. 19991 ). UKIRT K-band 



and HST-WFPC2 fiNiemela et al. 19981 ) optical imaging. Notable is the dominance of 
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the non-thermal component, interpreted as synchrotron radiation, between the thermal 
components of the stars, which points towar ds a colliding wind s cenario operating in this 



system. High-resolution X-ray observations (jPittard et al. 20021 ) revealed an extended 



emission region with its centroid located in agreement with the location of the bow shock 
between the stars. The binary separation at 5 GHz was determined to 643 ± 157 mas, 
corresponding to 417 AU/ cosi at a distance of 650 pc, where i is the inclination of the 
system, [i = 90° corresponds to the observer being located in the plane of the stars.) Here, 
the WR-star is separated from the non-thermal source by 575 ± 15 mas, while the B-star's 
projected distance to the synchrotron source is 68 mas. This is found i n good agreement with 



a wind-mom entum ratio of y 



(Mu 



ooJOB 



ooiWN 



0.011 ±0.016 dWiUiams et al. 1997h 



tor] = 0.02 JPittard et al. 2002h . From luminosit y (Lhni = 5 • lO^L^: 



Morris et al. 200 



Crowther et al. 1997) measurements the 



and spectral (effective temperature Teg = 28500 K: 
spectral type of the main sequence sta r has been estimate d to BO. 5 flWilliams et al. 19971 ). 
and more recently proposed to 05-71 (iLepine et al. 200 ll ) from HST STIS spectroscopy. 
Mass loss rates for the WN star has been derived from the VLA measurements by 
Churchwell et al. (1992). This was later improved by Morris et al. (2000), who found, by 
comparing his ISO images and optical data with non-LTE line models, wind filling factors 
/ of order 0.1. Because M oc Mtmc/y/f this leads to a lowering of th e WR-star loss rate 



to MwR = 2.4 ■ lO-^Mo yr"^ for a wind velocity of vwn = 950 km/s flMorris et al. 20001 ). 
Setia Gunawan (2001) then used a wind-momentum ratio (of ~ 0.014) to conclude on the 
mass-loss of the main sequence companian: Mqb = 4 ■ 10~^MQyr~^ for vqb = 800 km/s. On 
the basi s of the binary separ ation WR 147 has to be classified as an extremely long period 



binary (jVan der Hucht 20011 ). yet it's precise period is still unmeasured. A lower limit of 
1350 yrs was derived by Setia Gunawan (2001). No values for the eccentricity of WR 147 
are documented in the literature so far. 



Although attempts have been undertaken to estimate the inclination of this system. 
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the implied range of values for this quantity is broad: fitting the extended image 
of the non-thermal radio com ponent to a bow shock model gave valu es of 45° ± 15° 



(VLA: 



Contreras et al. 1999h . 30° ± 15° (VIA: 



Contreras et al. 20041) . 41° (MERLIN: 



WiUiams et al. 19971 ) and < 30° fMERLIN: lOoughertv et al. 20031 ). It was conjectured that 



the different antenna characteristics and the treatment of the system noise may influence 
the modeling procedure noticable. Practically, this means that the inclination of WR 147 
is at least highl y uncertain. Modeli ng the extended Chandra image required an inclination 



angle of ~ 60° (jPittard et al. 20021 ). Consequently, we will use this quantity in this work 
as a free parameter, in an attempt to constrain its value from observational limits on the 
non-thermal high energy component of WR 147. 

At present, only the radio band observations provided signatures of a non-thermal 
particle population in massive star systems, and data have been obtained in a number of 
measurements. We use the summary table from Setia Gunawan (2001) for the purpose 
of this work. After extracting the thermal component from the total measured one, one 
finds a power-law between 1 -10 GHz with spectral index a.^ ^ 0.43 — 0.5 {S^ oc i/""'"; 



Setia Gunawan et al. 2001a 



3) 



Dougherty et al. 20031 ). and with a low-energy turn-over at 



~ 0.5 — 1 GHz. A note of caution must be put here: the different instrument setups and 
capabilities, analysis methods, etc. lead to up 50% flux differences at the various frequencies. 
This applies in particular to the high-frequency data points at 22 GHz and 43 GHz, which 
additionally possess by themselves large measurement uncertainties. Suggestions (e.g, 



Dougherty et al. 20031 ) for a hi gh frequency turnover at these energies are there fore not to 



be considered established (e.g., 



Setia Gunawan et al. 2001a 



Pittard et al. 2006i ). 



For the present work we use the synchrotron energy spectrum S,y oc z/^°'^ with a cutoff 
well above 50 GHz. This corresponds to an E^"^ power law particle spectrum extending 
beyond a few 10^ MeV for mCauss magnetic field strengths. 
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Flux variations in the light curve of WR 147 have be en reported occasionally 



Churchwell et al. 19921) . however, could not be confirmed (ISetia Gunawan et al. 2001a 



Skinner et al. 19991 ) . Instead, any flux variability has rather been found to occur on long 



time scales with variability amplitudes not in excess of < 15% (jSetia Gunawan et al. 2001a 
possibly indicating wind clumping. 



3. The model for non-thermal high-energy radiation in binary systems of 

massive stars 

The present work utilizes the colliding wind model of Reimer et al. (2006). Here the 
supersonic winds of a massive main sequence star and its companion flow radially from 
each star until they collide to form a region of shocked, hot gas where their ram pressure 
balances, separated by a contact discontinuity with a forward and reverse shock. It is this 
shock region where particles may be accelerated out of the pool of thermal wind particles 
with acceleration rate E/E = a = Vg|(cj. — l)/(3crKa) (with /ta the corresponding diffusion 
coefficient, Vj^h the sho ck velocity and Cr the shock compression) by the Fermi mechanism 



(lEichler fc Usov 19931 ). Behind the shock the gas expands sideways from the wind collision 



region al ong the contact sur face with a convection velocity V being a fraction (set to 0.5 



here; see iReimer et al. 20061 ) of the wind velocity at the stagnation point. Both, diffusion 
and convection by the stellar winds operate in the collision region with diffusion being 
dominant close to the stagnation point (named "acceleration region" in the following), and 
convection determining the particles' transport further away from the downstream region 
(called the "convection region"). For simplicity the geometry of the collision region, being 
the site of the non-thermal radiation component, is adopted as cylinder-like. This is a 
sensible approximation specifically for high energy particles, the focus of this work, since 
most relativistic particles will have lost already a significant amount of their original energy 
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when entering the convection region. 

The location of the shock front with respect to the main sequence OB star, xqb, is 
determined by the balance of the ram pressure of both winds that have generally reached 
their terminal velocities (foB, "^^wn) in long-period binary systems: xqb = 1+^-^ where 
71 = ^OBVOB g^j^^ £) ^Yie separation of the binaries, Mor, Mwn the mass loss rates of the 
WN- and OB-star, respectively. 

For a given stellar surface magnetic field Bg the field strengt h B in the collision reg ion is 



assumed to follow the well established magnetic rotator theory 



Weber fc Davis 1967). and 



adopting a t ypical rotation velocity for early type stars of ~ O.lVnr, ( 



Conti fc Ebbets 1977 



1 r 

Pennv 19991. Estim ates for Bg range from below 100 G (iMathvs 19991 ) up to 10^ G 



Ignace et al. 19981 ). In Reimer et al. (2006) we used 5, = 30 G for WR 147 which leads 
to approximative equipartition field values at the shock location. For the present work we 
shall consider a sensible range for the values of Bg (see Sect. 4). 

The plane of the binary system is inclined by an angle i with respect to the observer 
{i = it/ 2 corresponds to an observer in the plane of the stars; see Fig. [1]), and (pB, the angle 
between the projected sight line and the line connecting the stars is a measure of the orbital 
phase $ of the system. Because the argument of periastron is not known for WR 147, we 
will define here periastron passage at $ = = 0b, and 0^ = for the WN-star being in 
front of the B-star along the sight line. In the co-rotating system centered on the OB-star 
the line-of-sight to the observer is described by the angles (pL with tan^^ = sin^^ cotz, 
and 9l with cos9l = coscpB sinz. 

The emitting electron spectra are calculated self-consistently by solving the steady 
state diffusion-loss equation and including wind expansion and all relevant radiative losses 
(synchrotron, IC and losses due to Coulomb scattering and relativistic bremsstrahlung) 
which balance the acceleration gains. IC losses often extend into the Klein-Nishina regime. 
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which requires a rigorous treatment of this dominating loss process. In consequence, 
Klein-Nishina loss dominated electron spectra cutoff at higher energies as compared to 
electron spectra that are treated with the Thomson approximation throughout the whole 
particle energy range. In the convection region the particles loose energy radiatively as 
well as through expansion losses while in the post-shock flow. This leads to a dilution 
of the particle density and to a deficit of high energy particles in this region. Because 
the stellar target photons for IC scattering arrive at the collision region from a prefered 
direc tion, the full an gular dependence of the cross section has to be taken into account 
(e.g., 



Reynolds 19821 ). The radiation yield of IC scattering depends on the scattering angle 
^ph = (angle between the directions of the incoming (from the OB-star) and outgoing 
photons), which obviously depends on the location of the scattering electron as well as the 
orbital phase. It is 

cos^ph = COS 9l cos 9 + sin^^ sin 9 cos(0 — 0^) . (1) 

This leads to anisotropy effects: the cutoff energy and flux variations of the IC luminosity 
are dependent on the line of sight into the wind as large scattering angles tend to produce 
more energetic photons than small scattering angles. A detailed description can be found 
in Reimer et al. (2006). 



4. Parameter study and their constraints 

The emission model of Reimer et al. (2006) for modeling a binary system of known 
geometry (orbital phase, eccentricity e, and period - or binary separation, the argument 
of periastron, and the inclination angle i with respect to the plane of the sky) possess five 
free parameters: main sequence star surface magnetic field, diffusion coefficient Kg, , shock 
compression factor Cr , convection velocity V and normalization of the injected particle 
population. These are contrasted by various physical and observational constraints: In order 
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to accelerate particles to the observed relativistic energies the diffusion coefficient must be 
low enough to allow acceleration to overcome the Coulomb losses. At the same time, the 
Bohm-limit for the diffusion coefficient shall not be violated. A further observational limit 
may be imposed by the observed maximum photon energy of any non-thermal component. 
Prom the highest observed frequency data point of the synchrotron component of WR 147, 
43 GHz, we deduce the maximum electron Lorentz factor 7e,max to be at least or greater 
than several 10^ for mGauss magnetic field strengths assuming roughly equipartition. To 
reach those energies diffusion must be at least as efficient as Kg, — 8- 10^^, 1.4 • 10^^, 2 • 10^^, 
4 • 102^ 1.4 • 10^^ and 14 • lO^^cm^ s"^ for system inclinations i = 0°, 30°, 45°, 60°, 75° and 
85° {Bg = 30G), respectively, and serve therefore as upper limits for These limits 
take into account the shift of the shock location away from the OB-star with growing 
binary separation for increasing inclination angles {D — i^proj/ cosi), and the corresponding 
dilution of the photon density there. The electron spectra would then cut off due to 
radiative losses. The highest electron energy is attained by decreasing the diffusion time 
scale up to the Bohm limit. This limit depends again on the system inclination i for an 
increasing binary separation and shock distance to the OB-star, and hence a decreasing field 
strength at the shock location, with rising i. The corresponding cutoff energies lie, within 
a factor 2-3, around several 10^ MeV {Bg — 30 G— 100 G). The narrowness of this range 
justifies the use of a fixed diffusion coefficient for all inclination angles without much loss of 
constraining power. In the following we therefore use = 3.2 • lO^^cm^ s^^ {Bg = 30G) or 
Ka — 0.9 • lO^^cm^ s~^ {Bg — lOOG), independent of the system inclination. We note that 
any value Bg > 10 G, together with the here assumed rotator theory, is in agreement with 
the lower limits for the field strength derived from the Razin-Tsytovich effect. 

The particle spectra normalization is limited by energy conservation and particle 
number conservation: the total energy in form of relativistic particles shall not be larger 
than the available total kinetic wind energy of the system, and the number of relativistic 
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particles injected into the system must not exceed the number of wind particles entering 
the acceleration zone. The final normalization is then fixed within these limits by the 
requirement to reproduce the observed synchrotron spectrum, and is therefore dependent 
on the chosen magnetic field strength. 

With these constraints and limits at hands we will in the following explore the effects 
of varying the system geometry upon the expected IC radiation at high energies. This 
will then allow us to compare our studied geometries with the observational data, and set 
constraints on the orbital parameters in the WR 147 system. 



4.1. Inclination effects on the particle spectra 

The excellent radio images gave a projected binary separation of ~ 417 AU for 
a distance of 650 pc. After de-projecting one obtains "true" star-star distances of 

= Z}pi.oj/(10^^ cosz) = 6. 25-71. 8cm for inclination angles i = 0° — 85°, respectively. 
For this range of i, and together with the suggested wind momentum ratio r] ~ 0.014 (see 
Sect. 2), the shock location xqb is then calculated to xqb/^O^^ cm= 0.66 — 7.60 where we 
find magnetic field strengths of 24.9 - 2.2 mG (5, = 30 G) and 59.1 - 5.1 mG (fi, = 100 G) 
for negligible field amplification (e.g. due to a run-away streaming instability in the 



acceleration region (IBell &: Lucek 200ll )). In Fig. [2] we show the importance of the various 
electron loss processes for a system inclination of i = 60° and Bg = 30 G as an example. IC 
losses determine a broad portion of the electron spectrum, until the synchrotron losses take 
over when Gompton losses are already deep in the Klein-Nishina regime. At low energies 
(less then a few tens of MeV), Goulomb losses dominate. The synchrotron spectra from 
WR 147 indicate an underlying electron spectrum with approximately E^"^. This spectral 



requirement is in agreement with the assumption of a strong s. 
velocity of ~ 1/2i;ob and anisotropic diffusion Ka = 'ii^d (see 



lock Cr = 4, a conv ection 
Reimer et al. 2006. for 
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details) which we shall use in the following. Fig. [3] shows the solution of the diffusion-loss 
equation in the acceleration region (of size tq ~ 8 x 10^'^ cm) for various inclination angles. 
Above ~ 10^ — 10^ MeV inverse Compton losses enter the Klein-Nishina regime. The 
particle spectrum finally cuts off due to the particle's gyroradius exceeding the size of 
the acceleration region tq rather than due to radiative losses. (The acceleration region is 
treated in our model as a leaky box, with a free escape boundary at tq. tq is the distance 
from the stagnation point where the diffusive time scale equals the convection and therefore 
independent of the magnetic field.) Since this limit scales with the magnetic field strength, 
the cutoff energy decreases with rising i as a response of the higher field density at the 
shock location. In Fig. H] we show an example of our computed synchrotron spectrum 
in comparison with the data. 0.36%-ll% (5, = 30G) or 0.27%-8.3% (B, = lOOG) for 
i = 0° — 90° of the kinetic wind energy of the OB-star, L„ ^ 8 ■ lO'^^erg s~^, must be used 
for electron acceleration to account for the observed synchrotron flux level. Here, the extent 
of the emission region was taken into account up to values of ~ 5 x 10^^ cm. 

4.2. Observing angle effects on the 7-ray spectra 

For anisotropic IC scattering, as is the case in the present systems, the resulting IC flux 
levels and maximum photon energies for a given electron Lorentz factor depend strongly 
on the scattering angle 6'ph. Here electrons see incoming photons from a variety of angles, 
which leads to raising or lowering the seed photon density accordingly in the electron rest 
frame after Lorentz transformation. Therefore the maximum scattered power goes into the 
6l = 180° direction. Similarily for the energy of the photons (in the Thomson regime): 
for small viewing angles the electrons are forward scattering the soft radiation that is 
less energetic in the electron rest frame than in the head-on case because of the Lorentz 
transformation. With Eq. 1 we flnd the maximum scattered power and photon energy when 
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9l = 180°, i.e. the WN star being behind the OB-star along the sight hne. At this angle the 
photons interact with the electrons primarily head-on, and end up therefore more energetic 
after the scattering event. Fig. [S] demonstrates the decreasing IC flux level and maximum 
scattered photon energy with declining viewing angle 9l between the sight line and the 
line connecting the stars for a fixed binary separation. From a given binary separation and 
the observed projected star-star distance one can infer the system inclination. Geometrical 
considerations [cosOl = cos^^sinz) then provide a range of possible viewing angles 6l- 
The viewing angle can reach higher values for large inclination systems. On the other 
hand, a large inclination increases the true binary separation, leading to a weaker magnetic 
(and diluted photon field) at the shock location, and therefore to lower maximum particle 
energies (see Figl3]) and Compton scattering rates. This interplay is demonstrated in 
Fig. the maximum gamma-ray energy decreases with increasing i as the maximum 
electron energy does (see Figl3]). For i = 85° the electron spectrum cuts off already at 
~ 3 ■ 10"^ MeV, too low to generate ~ 100 GeV photons in this setting. 

Our model-predicted fluxes are compared to typical sensitivity limits of instruments 
operating in the indicated energy range, in addition to the recently obtained INTEGRAL 
and MAGIC upper limits. For the binary separation that corresponds to i > 80° the 
MAGIC upper limit has no constraining power on the system geometry, irrespective of 
orbital phase and eccentricity. Furthermore, the overall flux level exhibits a declining trend 
with increasing i, diminishing also the constraints imposed by the INTEGRAL limits upon 
the model settings if the observer is located sufficiently close to the plane defined by the two 
stars. Specifically, we find the following constraints from the X-ray and 7-ray upper limits 
for WR 147 including the geometrical constraints: the MAGIC upper limit allows system 
settings with i > 80°, or for i ~ 70° — 75° the allowed range 61 increases sucessively from 
< 20° to < 50°, all assuming Bs = 30 G. The hard X-ray limits constrain the geometry to 
Ol < 40° for i < 70° and Ol < 50° - 60° for i = 75° - 85°. The angle $5 which measures 
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the orbital phase for a given system eccentricity (and hence orbital period) is determined by 
the knowledge of 9^ and i. The corresponding orbital phases (for an argument of periastron 
of zero) that meet the derived (0^, i) constraints range from to ~ 0.4 depending on the 
system eccentricity. Accordingly, the orbital period can be as low as ~ 10^ years or as high 
as ~ 10^ years. Adopting higher fields, e.g. Bg = 100 G, the range of allowed settings 
increases further: i > 70°, or for i around ~ 60° the allowed 6l ranges from < 30° — 40°. 
Also here the hard X-ray limits are only more constraining than the MAGIC limit for 
systems at large inclination angles: if i = 75° — 85°, the viewing angle must be constrained 
to 9l < 80° — 100°. The corresponding orbital phases are mostly close to 0, but can 
occasionally reach values up to ~ 0.5, and periods are again from ~ lO^years to ~ lO^years, 
depending on WR 147's eccentricity. Any estimate for the system eccentricity given the 
observational constraints set by the presented analysis is therefore coupled with the orbital 
phase. 

In summary, from our analysis it appears that large inclination angles are preferred 
for the WR 147 system if particle acceleration meets the exploited optimal conditions. On 
the other side, if WR 147 is barely inclined with respect to the plane of the sky, particle 
acceleration can not be as efficient as anticipated here. 

This study demonstrates the potential of measurements within the IC component 
in colliding wind binary systems who show clear signatures of synchrotron radiation, for 
constraining geometrical quantities which otherwise are often difficult to gather, specifically 
in long period systems. For geometrically well known systems, it offers the possibility to 
constrain the mostly rarely understood particle acceleration efficiency. Clearly, for WR 147 
the Fermi LAT energy range and expected sensitivitjo over a considerably long mission 
timeline seems promising to ultimately put rather decisive constraints on this system's 



""^http: / / www-glast.slac.stanford.edu/software/IS / glast_lat_performance.htm 



-17- 



geometry. In addition, due to its good sensitivity over several orders of magnitude of energy, 
Fermi LAT may probe the cutoff energy of the putative high-energy photon component, 
thereby providing further constraints on the particle acceleration efficiency in wind-wind 
collision regions. 

Long integration times texp in short period (T) binaries smear out any flux and spectral 
variations from a varying sight-line with orbital phase, resulting in a lack of the required 
observables (flux/spectrum for a given phase) to pin down geometrical system parameters. 
In long period binaries with T 3> texp on the other hand, significant orbital fiux variations 
during the exposure time are not expected. Therefore long-duration system geometries can 
be deduced in this case, even if long integration times are required. If in addition the orbital 
period, or orbital phase at the time of observation, is known, constraints on the system 
eccentricity can be deduced. Such observations performed at several phases and sufficient 
equally distributed throughout the system's orbit may intensify the constraining power of 
this method. 



4.3. 7-ray opacity 

Above several tens of GeV 77- pair production can a priori not be neglected in the 
dense UV environment of massive stars. Since the corresponding cross section is dependent 
on the interaction angle, the expected attenuation therefore relies on the system geometry, 
most important, on the distance of the shock location with respect to the main sequence 
star. For WR 147, being an extremely long-period binary with a presumably small 
eccentricity, the photospheric UV radiation field at the location of the collision region is 
found low enough (~ 10^~^ cm~^) to allow here the neglect of 7-ray absorption due to 
photon-photon coUisions. The value of r^^ < 0.015 is not exceeded in any of the considered 
cases, corresponding to a fiux throughput of at least 98% for WR 147. 
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5. Summary and conclusions 

Despite its close proximity the system geometry (including inclination, period, 
eccentricity) of the very-long period binary WR 147 is yet not well known. It belongs 
to a group of massive binary systems whose radio emission shows a non-thermal 
component. From radio imaging the location of this component has been associated 
with the wind collision region, and the projected binary separation is well constrained. 
WR 147 is one of the p rime candidates for gamma ray production in this collision zone 



(IReimer fc Reimer 20071 ). VHE observations with the MAGIC telescope have been made in 
2007 for 30.3 hours of good data, which resulted in upper limits at -F>snGev < 1.1 -10"^^ cm~ ^ 



s~\ i^>330GcV < 3.1 ■ 10-^2 cm-2 s~^ and F>6oogcV < 7.3 ■ 10^^^ cm^^ s~^ (!Aliu et al. 20081 ). 
Further observational constraints of WR 147's non-thermal component are coming from 
INTEGRAL upper limits at hard X-rays and the GHz frequency radio synchrotron 
measurements. In this work we explored the potential of gamma ray measurements to 
constrain the system geometry by means of an inverse Compton model for high energy 
photon production. For this purpose we assume a high efficiency of particle acceleration 
that generally allows the production of GeV photons. Within this scenario we find: 



• For the chosen rapidness of electron acceleration the particle spectra are cut off due 
to the size constraint of the acceleration region rather than due to radiative losses. In 
this case the cutoff energy scales with the magnetic field strength in the acceleration 
region, and we find the highest cutoff energy for low inclination systems. 

• The MAGIC upper limits at energies above 80 GeV therefore lacks constraining power 
for large inclinations > 80° for WR 147. 

• Provided particle acceleration acts sufficiently effective to allow the production of GeV 
photons through inverse Compton scattering, then our analysis indicates a preference 
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for WR 147 to possess a large inclination angle. However, if i is small in WR 147, 
electron acceleration can not be as efficient as anticipated here. 

• Increasing the magnetic field strength in the collision region widens the range of 
parameter values for the WR 147 system geometry. 

• The synchrotron data from WR 147 imply the existence of > 10^ MeV electrons, and 
therefore the existence of at least hard X-ray photons produced by those particles. 
Sensitive hard X-ray/soft gamma ray measurements have therefore the potential to 
constrain the WR 147 system geometry, independent of the assumed acceleration 
efficiency. 

• Gamma ray attenuation for this very-long period binary system is not anticipated 
to cause distinctive observational signatures in the high energy spectrum. Any 
phase- locked flux variations at 7-ray energies are therefore likely to be attributed to 
the anisotropic nature of the IC process. 

With a set of eight free parameters (system inclination and eccentricity, orbital phase 
and period, OB-star surface magnetic field, diffusion coefficient, shock compression factor 
and convection yelocity) within the framework of the IC model for high energy photon 



production (e.g. iReimer et al. 20061 ). to be confronted with the observationally constrained 
quantities (normalization and spectral index of synchrotron spectrum, projected binary 
separation, MAGIC and INTEGRAL upper limits) it is clear that an in-depth determination 
of the so far barely known geometry of the WR 147 system must await further information 
of from measurements of the non-thermal high-energy emission component. 



This work is supported by the National Aeronautics and Space Administration under 
contract NAS5-00147 with Stanford University. We thank Diego Torres for reading this 
manuscript and valuable discussions, and the referee for a constructive report. 



-20- 
REFERENCES 

Abbott, D.C., Bieging, J.H., Churchwell, E. & Torres, A.V. 1986, ApJ, 303, 239 
Aharonian, F., et al. (The HESS collaboration), 2005a, A&A, 442, 1 
Aharonian, F., et al. (The HESS collaboration), 2005b, Science, 309, 746 
Aharonian, F., et al. (The HEGRA collaboration), 2005c, A&A, 431, 197 
Aharonian, F., et al. (The HESS collaboration), 2007, A&A, 467, 1075 
Albert, J., et al. (The MAGIC collaboration), 2007, ApJ, 665, 51 
Albert, J., et al. (The MAGIG collaboration), 2008, ApJ, 675, L25 
Albert, J., et al. (The MAGIG collaboration), 2006, Science 312, 1771 
Aliu, E., et al. (The MAGIG collaboration), 2008, ApJ 685, L71 

Barlow M.J., 1982, in: "Wolf-Rayet stars: Observations, physics, evolution", Mexico, 
lAU-Symp. 99, 149 

Barker P., 1986, in: "Workshop on the Connection between Nonradial Pulsations and 
Stellar Winds in Massive Stars", Boulder, PASP, 98, 44 

Bednarek, W., 2005, MNRAS, 636, 46 

Bell, A.R. & Lucek, S.G., 2001, MNRAS, 321, 433 

Benagha, P. & Romero, G.E. 2003, A&A, 399, 1121 

Benagha, P., Romero, G.E., Stevens, I.R. & Torres, D.F. 2001a, A&A, 366, 605 
Benagha, P., Cappa, G.E. & Koribalski, B.S. 2001b, A&A, 372, 952 



Casse, M. & Paul, J.A. 1980, ApJ, 237, 236 

Cassinelli, J. P. 1979, Annual review of astronomy & astrophysics, 17, 275 
Chen, W. & White, R.L. 1991, ApJ, 366, 512 

Churchwell, E., Bieging, J.H., van der Hucht, K.A., et al. 1992, ApJ393, 329 
Conti, P.S. & Ebbets, D. 1977, ApJ, 213, 438 

Contreras, M.E., Rodriguez, L.F., Tapia, M. et al. 1997, ApJ, 488, L153 
Contreras, M.E. & Rodriguez, L.F., 1999, ApJ, 515, 762 

Contreras, M.E., Montes, G. & Wilkin, P.P., 2004, Revista Mexicana de Astronomia y 
Astrofisica, 40, 53 

Crowther, P.A. 1997, in: The Effective Temperature of Hot Stars, ed. T. R. Bedding, A. J. 
Booth, & J. Davies (Dordrecht: Kluwer Academic Pubhshers), 137 

De Becker, M., Rauw, C, Pittard, J.M. et al., 2007, A&A, 472, 905 

Donati, J.-F., Babel, J., Harries, T.J. et al. 2002, MNRAS, 333, 55 

Dougherty, S.M., WiUiams, P.M., van der Hucht, K.A. et al. 1996, MNRAS, 280, 963 

Dougherty, S.M., Williams, P.M. & Pallaco, D.L. 2000, MNRAS, 316, 143 

Dougherty, S.M., Pittard, J.M., Kasian, L. et al. 2003, A&A, 409, 217 

Dougherty, S.M., Beasley A.J., Claussen, M.J., et al. 2005, ApJ, 623, 447 

Eichler, D. & Usov, V. 1993, ApJ, 402, 271 

Funk, S., Reimer, O., Torres, D. F. & Hinton, J. A., 2008, ApJ, 679, 1299 



-22- 

Gayley, K.G., Owocki, S.P. & Cranmer S.R. 1997, ApJ, 475, 786 
Igance, R., Cassinelli, J.R & Bjorkman, J.E. 1998, ApJ, 505, 910 

Kaul, R.K. & Mitra, A.K., 1997, in: Proc. 4th Compton Symposium, Eds: C.D.Dermer, 
M.S.Strickman, J.D.Kurfess, AIR Conf. Rroc. 410, 1271 

Lamers, H.J.G.L.M. & Casinelli, J.R, "Introduction to Stellar Winds", Cambridge 
University Rress, 1999 

Lepine, S., Wallace, D., Shara, M.M., Moffat, A.F.J., et al., 2001, AJ, 122, 3407 

Maeder, A. & Meynet, G. 1987, A&A, 182, 243 

Malkov, M.A. & Volk, H.J. 1995, A&A, 300, 605 

Marchenko, S.V., Moffat, A.F.J., Ballereau, D., et al., 2003, ApJ, 596, 1295 

Mathys , G. 1999, in: "Variable and Non-spherical Stellar Winds in Luminous Hot Stars", 
eds. Wolf, B. et al.. Lecture Notes in Rhysics, 523, 95 

Moffat A.F.J., 1996, In: Vreux J.M., Detal A., Praipont-Caro D., Gosset E., Rauw G. 
(eds.) Wolf-Rayct Stars in the Framework of Stellar Evolution. Proc. 33rd Liege 
International CoUoq., Univ. Liege, Inst. d'Astrophys., Liege, p. 199 

Montmerle, T. 1979, ApJ, 231, 95 

Morris, P.W., van der Hucht, K.A., Croiwther, P.A., et al., 2000, A&A, 353, 624 

Miicke, A. & Pohl, M. 2002, In: Moffat A.F.J. & St-Louis, N. (eds.) "Interacting winds 
from massive stars" , ASP conf. series, 260, 355 

Niemela, V.S., Shara, M. M., Wallace, D.J., et al. 1998, AJ, 115, 2047 

Penny, L.R. 1996, ApJ, 463, 737 



-23- 

Pittard, J.M., Stevens, I.R., Williams, P.M. et al. 2002, A&A, 388, 335 

Pittard, J.M., Dougherty, S.M., R. F. Coker, R.F. et al., 2006, A&A, 446, 1001 

Pittard, J.M. & Dougherty, S.M., 2006, MNRAS, 372, 801 

Rauw, G., Blomme, R., Waldron, W.L. et al. 2002, A&A, 394, 993 

Reimer, A., Pohl, M. & Reimer, O., ApJ, 644, 1118 

Reimer, A., Proceedings of the workshop "Clumping in Hot Star Winds", eds. W.-R. 
Hamann, A. Feldmeier and L. Oskinova, pp.187 

Reimer, A., Reimer, O. & Pohl, M., 2007, Ap&SS, 309, 351 

Reimer, A. & Reimer, O., 2007, AIP Conf. Proc, 921, 217 

Reynolds, S.P. 1982, ApJ, 256, 13 

Romero, G.E., Benagha, P., Torres, D.F. 1999, A&A, 348, 868 
Schild, H., Gdel, M., Mewe, R., et al. 2004, A&A, 422, 177 

Setia Gunawan, D.Y.A., de Bruyn, A.G., van der Hucht, K.A., et al. 2001, A&A, 368, 484 

Setia Gunawan, D.Y A., van der Hucht, K.A., WiUiams, P.M. et al. 2001, A&A, 376, 460 

Skinner, S.L., Itoh, M., Nagase, F., et al., 1999, ApJ, 524, 394 

Skinner, S.L., Zhekov, S.A., Giidel, M. et al, 2007, MNRAS, 378, 1491 

Stevens, I.R., Blondin, J.M. & Pollock A.M.T., 1992, ApJ, 386, 265 

Stevens, I.R. & Pollock, A.M.T., 1994, MNRAS, 269, 226 

Usov, V.V. 1992, ApJ, 389, 635 



Usov, V.V. & Melrose, D.B., 1992, ApJ, 395, 575 

Van der Hucht, K.A., 2001, New Astronomy Reviews, 45, 135 

Weber, E.J. & Davis, L.Jr., 1967, ApJ, 148, 217 

White, R.L. 1985, ApJ, 289, 698 

White, R.L. & Chen W. 1992, ApJ, 387, L81 

White, R.L. & Becker, R.H. 1995, ApJ, 451, 352 

WiUiams RM., van der Hucht, K.A., Pollock, A.M.T., et al. 1990, MNRAS, 243, 662 
WiUiams P.M., Dougherty, S.M., Davis, R.J., et al, 1997, MNRAS, 289, 10 
Woosley, S. 1993, ApJ, 405, 273 
Zhekov, S.A., 2007, MNRAS, 382, 886 



This manuscript was prepared with the AAS macros v5.2. 



- 25 - 




Fig. 1. — The binary orbit of WR 147 for an assumed inclination of 60° and eccentricity 
e = 0.7. For details see text. 
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Fig. 2. — Energy loss rates in the colliding wind region of the WR 147 binary system, assumed 
to be inchned by 60° with respect of the sky: IC scattering (thin sohd lines) in the Thomson 
regime (TL) and Klein-Nishina regime (KN), synchrotron radiation (dashed-triple-dotted 
line), relativistic electron-ion bremsstrahlung (dashed-dotted line) and Coulomb interactions 
(dashed line) in comparison to the acceleration rate (thick solid line) in the acceleration 
zone. The dashed line represents the IC loss approximation used in this work. Parameters 
are: Lboi,oB = 5 • IO^Lq, ct = 6.6 eV, Mqb = 4 • lO-^Moyr-^, Mwn = 2.5 • IQ-^M^yr-^, 
vqb = 800km/s, -^wn = 950 km/s, B, = 30G, D = 417AU/cosi ^ 12.5-10i5cm, xqb ~ O.ID, 
B ^ 12mG, nph,T ~ 5.5 • 10W-^ Nh ~ 3.4 • lO^cm'^ k„ = 3.2 • lO^^cm^s-^ V ^ 399 km/s. 
To ~ 2 • 10^ s, ro 8 • lO^^cm. 
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Fig. 3. — (Non-normalized) steady-state electron spectrum for WR 147 for an inclination 
angle of i = 0°, 30°, 45°, 60°, 75° and 85° (upper to lower curves, each artificially separated 
by one order of magnitude), Qo — 1 and assuming a surface magnetic field for the B-star of 
30 G. See text for further parameters and discussion. 
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Fig. 4. — Synchrotron spectra for WR 147 from the i = 60° electron spectrum as shown 
in Figj3] and using the ^-approximation in comparison to the observed non-thermal radio 
emission taken from Setia Gunawan (2001a). The bulk of the synchrotron component is 
hidden below the photospheric UV radiation from the B-star. The thermal X-ray emission 
from the hot shocked gas in the collision region is not shown here. 
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Fig. 5. — IC spectra of WR 147 for a stellar separation of 6.25 • lO^^cm and varying angle 
9l = 0°, 10°, 20°, 30°, 40°, 50°, 60°, 90°, 120°, 150°, and 180° (lower to upper curves) to 
the observer. 0.36% of the B-star (with an assumed surface magnetic field Bg — 30G) wind 
kinetic energy has been injected into the acceleration region in form of electrons. 7-ray 
absorption by 77 pair production is found negligible at all energies. 
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Fig. 6. — IC spectra of WR 147 for an inclination angle of i = 30° (corresponding to a 
stellar separation of 7.22 ■ lO^^cm) and geometrically possible angles 6l = 60°, 90° and 120° 
(lower to upp er curves) to the observer and compared to Fermi-LAT's sensitivity of the 
inner Galaxy (IFunk et al. 20081 ). 0.36% of the B-star (with an assumed surface magnetic 
field Bs = 30G) wind kinetic energy has been injected into the acceleration region in form 
of electrons. 7-ray absorption by 77 pair production is found negligible at all energies. The 
MAGIC (purple arrows) and INTEGRAL (red arrows) upper limits rule out this inclination 
angle for WR 147 for the assumed free parameter set (see text for details), or the electron 
spectrum at the highest energies falls off more steeply than anticipated, possibly due to 
particle acceleration being not as efficient as anticipated in this work. 
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Fig. 7. — Same as Fig. [6] but for an inclination angle of i = 45° (corresponding to a stellar 
separation of 8.85 • lO^^cm) and geometrically possible angles 9l = 50°, 60°, 90°, 120° and 
130° (lower to upper curves). 0.55% of the B-star wind kinetic energy has been injected into 
the acceleration region in form of electrons. 7-ray absorption by 77 pair production is found 
negligible at all energies. The MAGIC (purple arrows) and INTEGRAL (red arrows) upper 
limits rule out this inclination angle for WR 147 for the assumed free parameter set (see text 
for details), or particle acceleration is not as efficient as anticipated in this work. 
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Fig. 8. — Same as Fig. [6] but for an inclination angle of i = 75° (corresponding to a stellar 
separation of 2.42-lO^^cm) and geometrically possible angles Ol = 20°, 30°, 40°, 50°, 60°, 90°, 
120°, 150°, and 160° (lower to upper curves). 2.23% of the B-star wind kinetic energy has 
been injected into the acceleration region in form of electrons. 7-ray absorption by 77 pair 
production is found negligible at all energies. The MAGIC (purple arrows) and INTEGRAL 
(red arrows) upper limits constrain the available parameter space to 6l < 50°. 
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Fig. 9. — Same as Fig. [6] but for an inclination angle of i = 85° (corresponding to a stellar 
separation of 7.18 ■ lO^^cm) and geometrically possible angles 6l = 10°, 20°, 30°, 40°, 50°, 
60°, 90°, 120°, 150°, and 170° (lower to upper curves). 11% of the B-star wind kinetic energy 
has been injected into the acceleration region in form of electrons. 7-ray absorption by 77 
pair production is found negligible at all energies. The INTEGRAL (red arrows) upper 
limits constrain the available parameter space to 9^ < 60°, the MAGIC upper limits does 
not provide any constraint. 



